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Abstract 

Canopy gaps created by wind-throw events, or blowdowns, create a complex mosaic of forest patches varying in 
disturbance intensity and recovery in the Central Amazon. Using field and remote sensing data, we investigated the short- 
term (four-year) effects of large (>2000 m ) blowdown gaps created during a single storm event in January 2005 near 
Manaus, Brazil, to study (i) how forest structure and composition vary with disturbance gradients and (ii) whether tree 
diversity is promoted by niche differentiation related to wind-throw events at the landscape scale. In the forest area affected 
by the blowdown, tree mortality ranged from 0 to 70%, and was highest on plateaus and slopes. Less impacted areas in the 
region affected by the blowdown had overlapping characteristics with a nearby unaffected forest in tree density (583±46 
trees ha -1 ) (mean±99% Confidence Interval) and basal area (26.7±2.4 m 2 ha -1 ). Highly impacted areas had tree density 
and basal area as low as 120 trees ha -1 and 14.9 m 2 ha -1 , respectively. In general, these structural measures correlated 
negatively with an index of tree mortality intensity derived from satellite imagery. Four years after the blowdown event, 
differences in size-distribution, fraction of resprouters, floristic composition and species diversity still correlated with 
disturbance measures such as tree mortality and gap size. Our results suggest that the gradients of wind disturbance 
intensity encompassed in large blowdown gaps (>2000 m 2 ) promote tree diversity. Specialists for particular disturbance 
intensities existed along the entire gradient. The existence of species or genera taking an intermediate position between 
undisturbed and gap specialists led to a peak of rarefied richness and diversity at intermediate disturbance levels. A diverse 
set of species differing widely in requirements and recruitment strategies forms the initial post-disturbance cohort, thus 
lending a high resilience towards wind disturbances at the community level. 
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Introduction 

Natural disturbances varying in size are suggested as major 
driver of tree species substitution across space and time, thus 
influencing vegetation structure, composition and diversity in 
forest ecosystems [1-4]. Disturbance dynamics in neotropical 
forests are considered to be dominated by small canopy gaps 
(<2000 m 2 ) from small treefall events resulting from small-scale 
abiotic disturbances [5-7] and biotic interactions [8]. In contrast, 
widespread tree mortality caused by large-scale exogenous 
disturbances such as forest blowdowns are often regarded as rare 
events. 

However, remote sensing studies confirm a broad range in the 
size and intensity of wind disturbances representing a continuum 
without simple distinctions between frequent and episodic events 



[9-11]. In addition, widespread tree mortality associated with 
blowdowns is more prevalent in the Central and Western Amazon 
than the Eastern Amazon [12-14]. This has revived a classical 
question [1,5,15,16] on the ecological importance of large canopy 
gaps (2000 m 2 ) and how these influence forest dynamics and tree 
diversity patterns at the landscape scale [11,14,17-20]. 

In Central Amazon forests, there is in fact a large gradient of 
gap sizes created by blowdowns, which leads to a spatially complex 
mosaic of successional pathways [10,20,21] and gap size distribu- 
tion varies locally, for example with topography. Studies outside 
the Amazon report that gap size distributions resulting from wind 
disturbance vary predictably with elevation [22-24], with a greater 
fraction of large contiguous blowdowns on exposed ridges and 
plateaus. Blowdowns can raze thousands of trees locally, and with 
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this create a few large gaps ( 2000 m 2 ) as well as many smaller tree- 
fall gaps (<2000 m 2 ) [9,10,25]. 

Large blowdowns and small tree fall gaps may differ substan- 
tially in environmental conditions and thus may represent very 
different starting conditions for tree regeneration. In single-treefall 
gaps, soil and organic layer disturbance due to uprooting and 
snapping trees, as well as changes in nutrients and water 
availability can influence both mortality and recruitment rates as 
well as other factors vital for the maintenance of tree populations 
and communities [6,15,26-28]. Large blowdowns can include a 
wide range of variability in disturbance severity [22,29], which 
may affect the structure and composition in the recovering forest 
[30-33]. 

As a first order effect, gap formation by adult tree mortality is 
considered a key process as it increases the local availability of 
resources, most notably of light [5,34]. With the resulting increase 
in high-light microsites [26] forest patches gain the potential to 
promote light-demanding (pioneers) tree species specialized in 
colonizing and occupying gaps by efficient dispersal and rapid 
growth [2,34—36]. Nonetheless, the same traits that convey high 
growth potential for such light-demanding or pioneer species (e.g. 
high enzyme activity, low self-shading, and low wood density and 
construction costs) weaken their competitive strength under 
undisturbed conditions. Such traits may prevent them from 
enduring low light conditions, casting shade and outcompeting 
shade-tolerant and slower-growing (late-successional) trees in the 
longer term, resulting in relatively short life spans [37-41]. Light- 
demanding species thus depend on frequent gap formation to 
escape local extinction [2,5,6,34]. However, this binary mosaic 
view - pioneer versus late-successional species in gaps versus old 
growth forest - may be too simplistic for understanding the 
influence of a complex gradient ranging from single tree fall gaps 
to large blowdowns on landscape-level dynamics of highly-diverse 
tropical rain forests. Given the complexity of the disturbance 
mosaic, tropical tree species may not fall into two distinct 
successional groups but rather form a continuum from pioneer 
to late-successional strategies [37,42-44]. For example, small gaps 
(<2000 m 2 ) opened up by mortality events of one or a few trees 
maybe too dark for light-demanding species and colonization may 
depend on alternative mechanisms, such as growth-release of the 
understory sapling bank (advanced regeneration), resprouting of 
damaged trees, and lateral expansion of surviving trees. All three 
mechanisms may also be relevant in larger gaps (>2000 m 2 ) 
especially given distance-related limitation of seed dispersion. 
[31,45-48], but here high-light conditions also strongly promote 
the emergence of pioneer species [33,49]. However, the threshold 
disturbance intensity below which colonization by light-demand- 
ing species fails is not known. Moreover, it is not known whether 
there are specialists for intermediate levels of wind disturbance. 

By combining a landscape-level approach with a detailed plot- 
scale analysis of high taxonomic resolution, we analyzed recruit- 
ment patterns along a large gradient of disturbance intensity and 
its potential relevance for maintaining tree species diversity in 
tropical forests. We hypothesize that the vast species pool of 
Amazonian forests contains species with strategies intermediate 
between the extremes of light-demanding and shade-tolerant 
species, thus enabling the forest vegetation to fully exploit the 
entire gradient of gap sizes and resource availability. We further 
hypothesize that a co-occurrence of light-demanding, shade- 
tolerant and intermediate species leads to a peak in diversity at 
intermediate disturbance intensities and gap sizes in early 
regeneration. To address these hypotheses we answered the 
following questions: 1. Are there predictable patterns of disturbance 
intensity distributions for large blowdowns, and do these patterns 



differ between topography classes (valleys, slopes, and plateaus)? 2. 
Do severe blowdowns exert/promote selective mortality effects at 
individual and species dimensions'? 3. How do pre-blowdown 
conditions, early vegetation responses and the relative importance of 
the regeneration mode (establishment from seed, resprouting and 
growth of survivors) interact with disturbance intensity? 4. What 
are the implications of these effects on community composition and 
species coexistence at the landscape scale? 

Materials and Methods 

Study sites 

We conducted the research at the Estacao Experimental de 
Silvicultura Tropical (EEST) (2°61'S, 60°20'W) of the Instituto 
Nacional de Pesquisas da Amazonia (INPA) and adjacent- 
contiguous area (ZF2), which is administered by the Super- 
intendencia da Zona Franca de Manaus (SUFRAMA) (2°56'S, 
60°26'W), Amazonas, Brazil (Figure la). In January of 2005 
storms propagating from southwest to northeast of Brazil caused 
large forest blowdowns across the central Amazon, including 
—2500 ha of forest in the Manaus region [9]. At the ZF2, forest 
not known to have been previously disturbed near the Rio 
Cuieiras (a tributary to the Rio Negro) was heavily disturbed, as 
identified by field surveys and Landsat images (Figure 1). 

The entire area is covered by terra firme forest. Prior to the 
2005 disturbance there was no evidence of human intervention for 
at least the previous 40 years. In our Landsat chronosequence 
there was no signal of large natural disturbances affecting the area 
since 1979. Mean monthly temperature in this region is 26°C with 
litde seasonal variation [50,51] and annual rainfall is about 
2600 mm, with a distinct dry season between July and September 
[52,53]. The local topography is undulating with a maximum 
altitudinal difference of about 140 m (40-180 m a.s.l). Upland 
plateaus with high clay content (Oxisols) are cut by slopes and 
valleys dominated by soils with high sand content (Spodosols) and 
subject to seasonal flooding. The drainage network flows to 
tributaries of the Rio Negro or directly to the Rio Cuieiras. The 
terra-firme forest is characterized by a closed canopy with high 
tree species diversity [54—57] and a dense understory with 
abundant acaulescent palm species in plateaus and canopy palm 
species in valleys [57,58]. In this region, forest compositional and 
structural variations are correlated with water stress, soil and 
topography [54,59-61]. 

Forest inventory 

We sampled three different patches direcdy affected by the 2005 
blowdown (hereafter referred to as disturbed forest). To assess the 
entire disturbance gradient we installed three pairs of transects, 
measuring 200, 600 and 1000 m length by 10 m width 
respectively, for a total of six transects (Figure lb - blue lines). 
We established 144 sub-plots of 10x25 m (total of 3.6 ha) within 
the 6 transects. At each sub-plot, we determined several measures 
of forest structure for trees with diameter at breast height (DBH) 
S10 cm. Diameters of live uprooted and damaged-snapped trees 
were usually measured above the DBH (1.3 m height) to avoid 
lesions and mechanical obstructions, such as trunks, branches, soil, 
etc. For species with buttressed and aerial roots, we measure the 
diameter just above these obstructions. 

Considering previous studies [26,46,62] and field observations, 
we identified resprouting stimulated by mechanical injuries in 
individuals following uprooting and snapping (partial or total 
rupture of the crown). For botanical identification we collected 
samples from at least one individual of each species and for 
subsequent re-measurement we tagged all trees. We carried out 
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Figure 1. Study areas at the confluence of the Rio Solimdes and Rio Negro, Amazonas, Brazil. Legend: (a) Landsat RGB composition of 
the studied areas (red inset); (b) sampled areas [short-wave infrared reflectance (red channel) indicate the 2005 blowdown tree mortality, measured 
by quantifying the differences on the no-photosynthetic vegetation (ANPV) fraction - SMA analysis]; (c) SRTM elevation model and topographic 
variation. 

doi:1 0.1 371 /journal.pone.01 0371 l.gOOl 



identification to the species level [63-65] when possible, by 
comparing the collected material with specimens at the INPA and 
IFAM (Institute Federal de Educacao, Ciencia e Tecnologia do 
Amazonas) herbariums, and also with the EEST botanical 
collection. Reproductive material (flowers and/or fruits) was 
collected when available, and added to the IFAM herbarium 
and to the EEST botanical collection. Our research did not 
involve endangered or protected species and no special permits 
were required to access the ZF2. We carried the forest inventory 
between June and September of 2009. Botanical samples were 
collected in two campaigns, in 2009 and 201 1. 

We used measurements from an undisturbed contiguous area 
(hereafter referred to as undisturbed forest) to compare with our 
data for the disturbed forest (Figure lb - green lines). In this area, 
two transect plots measuring 20x2.500 m (10 ha) were installed in 
1996 [66] and forest structure and dynamics are monitored in 250 
sub-plots (20x20 m) using measures of growth, recruitment and 
mortality made every two years since 2000 [52,67]. We used data 
from 196 (total 7.8 ha) sub-plots in an area not affected by clouds 
in both 2004 and 2005 Landsat images, and estimated structural 
parameters at these sites in 2008. We used census data from these 
sites for the years 2004-2006 to obtain background forest 
mortality rates for an area not directly affected by the 2005 
blowdown (as confirmed using our Landsat ANPV maps) [10]. 
Although much of the Amazon experienced drought in 2005 with 
potential effects on mortality [68,69], the region we studied was 
not unusually dry during this period [9,53]. 

The transects established in disturbed and undisturbed forests 
were perpendicularly oriented in north-south and east-west 
directions (Figure lb). These transects cross several toposequences 
in order to represent the floristic and structural variation between 
plateaus, slopes and valleys (Figure lc). For classifying sub-plots 
into topographic classes we used field observations, apparent soil 
texture, slope and altitude [56,57,66]. Slope was measured with a 
clinometer and sub-plot altitude data was assessed by using a 
Global Positioning System (GPS) receiver (Trimble Nomad 900) 
and a digital elevation model (Shutde Radar Topographic Mission 
-SRTM) (http://glovis.usgs.gov). Higher altitude areas (>80 m 
a.s.l.) with lower slope (~0— 10°) and higher clay content soils were 



classified as plateaus. Areas with slope >10°, mixed sand/clay 
texture soils, and located in interfluve areas were classified as 
slopes. Areas with low elevation and relief with sandy soil and 
direct contact with perennial or intermittent water streams were 
classified as valleys (Figure 1). In the Manaus region, valleys can 
be also characterized by the high abundance of some palm species 
(e.g. Mauritiella aculeata (Kunth) Burret, Oenocarpus bataua 
Mart, and Manicaria saccifera Gaertner) and for being partially 
inundated during the rainy season. In the disturbed forest we 
classified 44 (1.1 ha), 73 (1.8 ha) and 27 (0.7 ha) sub-plots in 
plateaus, slopes and valleys, while in the undisturbed forest 
plateaus, slopes and valleys represented 49 (2.0 ha), 96 (3.8 ha) and 
51 (2.0 ha) sub-plots, respectively. 

In both disturbed and undisturbed forests we counted dead trees 
and measured their DBH. In the disturbed forest, we classified 
mode of tree death related to the blowdown as snapped, uprooted 
or standing dead according to previous studies [9,25,62]. For both 
disturbed and undisturbed forests, we estimated tree density (trees 
ha ) and basal area (m 2 ha ). For the disturbed forest, we 
additionally estimated mean DBH (cm) and wood density (g cm 3 ) 
for each sub-plot. We compiled species wood density from 
available literature [70-73] to check for possible correlations 
between disturbance intensity measures and wood density 
variation. For species for which data were not available or the 
identification was only possible at the genus level, we used genus- 
level obtained by calculating the mean wood density of local 
congeneric species. For six trees (snapped or without leaves), we 
used family-level obtained by calculating the mean wood density of 
the recorded species from die respective family. 

Satellite data and disturbance parameters 

We used Landsat-5 thematic mapper images (Path 231, Row 
062), obtained from Brazil's National Institute for Space Research 
(INPE, http://www.dgi.inpe.br/CDSR/) and the United States 
Geological Survey (USGS, http://glovis.usgs.gov) to assess distur- 
bance intensity over the disturbed area, directly affected by the 
2005 blowdowns. We first georeferenced (400 control points per 
image) the images using the NASA Geocover data set as a 
basemap (https:/zulu.ssc.nasa.gov/mrsid/) and several reference 
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points across the EEST route were collected using a GPS (Trimble 
Nomad 900); we applied a mask to remove clouds, cloud shadows, 
land use, and water bodies; and removed smoke and haze 
contamination with the Carlotto technique [74]. Further, we 
carried out a spectral mixture analysis (SMA) [75] to determine 
the fractions of green vegetation (GV), nonphotosynthetic vege- 
tation (NFV) and shade. We normalized the pixels without shade 
as GV/(GV + NPV) and NPV/(GV + NPV), resulting in shade- 
normalized GV and NPV fractions [76]. The ANPV reflects 
changes in nonphotosynthetic vegetation (NPV, wood, dead 
vegetation and surface litter) before and after the storm and was 
calculated as NPV2005-NPV2004. We carried out SMA analysis 
with the Environment for Visualizing Images software (ENVI, 
ITT Industries, Inc, Boulder CO, USA). 

We assessed disturbance intensity in the disturbed forest by 
estimating three disturbance parameters on a per-pixel basis 
(30x30 m): sub-plot mortality, mortality in the area surrounding 
the sub-plot (hereafter referred to as neighboring mortality) and 
gap size; these three disturbance parameters derived from tree 
mortality measures (in percentage) obtained from a mortality 
estimation model which has the ANPV as predictor 
[Y= 103.22*ANPV-3.32] (^ = 0.8 and P<0.001) [9]. We used 
the four corner coordinates (northeast, northwest, southeast and 
southwest) to rasterize polygons representing our sub-plots on the 
ANPV image. The sub-plot's ANPV was obtained by calculating 
the area weighted mean ANPV value of those pixels included in its 
respective rasterized polygon. To increase the accuracy of this 
estimation, we resampled the ANPV image from its native 
resolution of 30x30 m to a 3x3 m grid size, which allowed 
improved estimation of the relative area within each pixel 
contributing to the total sub-plot area. 

We applied the sub-plot's weighted mean ANPV in our 
mortality estimation model [9] to assess tree mortality at sub- 
plot level in the two forests. For the disturbed forest we calculated 
the neighboring mortality with the same equation by applying the 
mean ANPV value of the eight pixels directly adjacent to that pixel 
which included the geographic coordinates of the sub-plot's central 
point. Additionally, we estimated gap size as the area of virtual 
polygons formed by groups of continguous pixels with ANPV& 
0.16. In these forests a ANPV^O.16 corresponds to a mortality 
rate of ~13% (~6 trees/pixel), which is conservatively above the 
minimum rate observed from smaller treefall gaps in local old 
growth forests (~2%) [67,77], and above the single-pixel mortality 
gaps obtained from field-measured mortality (~5%) [10,25]. In 
our mortality model [Y = 103.22*ANPV-3.32], ANPV<0.032 
produced negative mortality estimates, which occurred exclusively 
in the undisturbed forest and valleys of the disturbed forest. 
Because the model was adjusted only for disturbed areas, negative 
ANPV values cannot be interpreted as real gain in PV (or forest 
regrowth). In order to avoid negative mortality estimates, we 
applied a threshold (ANPV<0.032) below which we assumed 
mortality was zero for sub-plots. Neighboring mortality and gap 
size were only estimated for the disturbed forest. We computed 
sub-plot and neighboring mortality by using tools from the raster 
[78] and maptooh [79] packages implemented in the R (version 
3.0.1) software platform [80]. We computed gap size with the 
SAGA software (version 2.0.8 http://www.saga-gis.org). 

Statistical analysis 

Disturbance patterns. To check for recognizable relation- 
ships between disturbance patterns and topography, we compared 
structural (tree density and basal area) and disturbance intensity 
(sub-plot mortality) of the different forests. For the disturbed forest, 
we tested for the effects of disturbance intensity measures (sub-plot 



mortality, neighboring mortality and gap size) on subsequent 
structural and floristic variation in the different topographic classes 
of the disturbed forest. 

To assess structural differences among topographic classes 
within the two forests, we classified sub-plots from the disturbed 
forest in two levels of disturbance. Non- and low-disturbance sub- 
plots (ANPV<0.16; hereafter referred as low disturbance) [10] 
were compared with sub-plots that experienced disturbance above 
this threshold (ANPV>0.16; hereafter referred to as high 
disturbance). We compared ANPV and structural measures of 
the two forests with factorial ANOVA. We used a Post-hoc Tukey 
test to evaluate significant differences among topographic classes. 
In disturbed forest, we related mode of death (snapped, uprooted 
and standing dead) among topographic classes with linear 
regressions and assessed whether mode of death is still recogniz- 
able and related to sub-plot mortality. With one-way ANOVA, we 
compared the mean DBH of dead trees among modes of death. 
Additionally, we tested for the effect of the disturbance intensity on 
subsequent structural variation by relating tree density and basal 
area to sub-plot mortality with linear regressions. The same 
analysis was applied to test for possible correlations between mean 
DBH of recorded dead trees and sub-plot mortality. 

Selective mortality patterns. To address our questions 
related to selective mortality effects at individual and species level, 
we compared the size distribution of trees from the disturbed with 
the undisturbed forest. For the disturbed forest, we described the 
variation of genera's importance at low and high disturbance 
levels, and related species diversity measures to sub-plot mortality. 

Due to the lower disturbance intensity that has been observed in 
valleys, and also to pre-blowdown distinct structural and floristic 
characteristics of the vegetation in the valleys, we only considered 
sub-plots established on plateaus and slopes for these analyses. To 
assess the influence of disturbance intensity on size distribution of 
trees and on the importance of genera, we also considered the 
classification described in section 2.4.1. We assessed the distur- 
bance intensity effects on the size distribution of trees using Chi- 
squared tests comparing the diameter distribution of live and dead 
trees within the disturbed and undisturbed forests, and also within 
the levels of disturbance described above. For the disturbed forest, 
we additionally related the mean DBH of live trees to sub-plot 
mortality. We described community demographic patterns at the 
genus level for low and high disturbance sub-plots by estimating 
the Importance Value Index (IVI), calculated as the sum of relative 
density, frequency and dominance (basal area) of congeneric 
species [81]. To assess the effects of disturbance intensity on 
diversity patterns, we related with linear and polynomial 
regressions, species richness, Shannon diversity and species 
rarefaction to sub-plot mortality. These three measures were 
calculated for groups of five sub-plots, ranked and grouped by 
their sub-plot mortality value (23 in total). With this we avoided 
the tree density effect observed in those sub-plots with high 
mortality values, and thus few trees (e.g. with three trees); and were 
able to look for patterns at the community level. 

Disturbance gradient effects and vegetation 
responses. To address our last two questions, on vegetation 
responses along the disturbance gradient and changes on 
community composition, we assessed the importance of resprou- 
ters and fast growing or pioneer species, and checked for variations 
in species composition with disturbance intensity measures. 

With linear regressions we assessed the effect of sub-plot 
mortality, neighboring mortality and gap size on the fraction of 
resprouters; and effect of sub-plot mortality on mean wood 
density. We tested with one-way ANOVA for wood density and 
size differences between resprouters and non-resprouters. Further, 
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we related the fraction of pioneer species basal area to sub-plot 
mortality and mean wood density; and related the fraction of 
pioneer species with neighboring mortality and gap size. We 
consider as pioneers, light-demanding/ fast-growing species typical 
for the region from the genera Cecropia, Conceveiba, Inga, Laetia, 
Miconia, Pourouma, Tachigali and Vismia. To reduce the 
influence of possible pioneer trees that were established before 
the blowdown, we only accounted for trees in these genera with 
DBH<25 cm. 

To test for the effect of disturbance intensity on floristic 
composition, we used nonmetric multidimensional scaling 
(NMDS) to reduce the dimensionality of the community species 
composition [82,83]. Ordinations were computed with two axes 
from one dissimilarity matrix generated by the abundance of all 
recorded species. We used Mantel-tests (P<0.001) to test for the 
significance of the NMDS analysis. Finally, sub-plot mortality, 
neighboring mortality and gap size were used as predictors in 
linear regressions to assess the effects of disturbance intensity 
measures on the floristic composition variation detected by the first 
NMDS axis. Finally, to look in detail at possible variations on 
genera composition driven by the disturbance intensity, we fit a 
cubic smoothing spline function (degrees of freedom = 3) relating 
the abundance of the 25 most important genera (based on IVI 
values) to sub-plot mortality. As in our previous analyses (see 
section 2.4.2) only sub-plots established on plateaus and slopes 
were used for these last two analyses. All statistical analyses were 
performed in R (version 3.0.1) software platform [80]. We used 
the vegan package [84] for estimating diversity indices, species 
curves and to process the NMDS analysis. All other tests were 
custom written. 

Results 

Disturbance patterns 

In the disturbed forest we sampled 1944 live and 363 dead trees. 
Our structural measures confirm a correlation between ANPV and 
tree mortality, and indicate that vegetation damage was partially 
controlled by topography. 

Tree density and basal area in the sub-plots (250 m 2 ) of the 
disturbed forest ranged from 3 to 33 trees and from 0.08 to 2.23 
m 2 , respectively. Mean tree density and basal area in the 
undisturbed forest was higher than in the low and high disturbance 
sub-plots (Table 1). Mean tree density and basal area in slopes 
from the high disturbance was lower than in slopes from the low 
disturbance and undisturbed forest (Table 1). 

In the disturbed forest, ANPV ranged from -0.213 to 0.696 
with an overall mean value of 0.176±0.037 (mean±99% CI), 
higher than from the undisturbed forest (0.038 ±0.0 16; ranging 
from —0.168 to 0.202). Estimated sub-plot mortality in the 
disturbed forest ranged from 0 to 69.9% with an overall mean 
value of 16.2±3.4%, higher than from the undisturbed forest 
(3.7±0.8%; ranging from 0 to 30.5). In the undisturbed forest, 
mortality rate computed from field data (2004-2006) was 
1.6 ±0.7% year and correlated positively but weakly with the 
estimated sub-plot mortality (r 2 = 0.04; P = 0.004; Pearson's 
r = 0.20). Plateaus, slopes and valleys from high disturbance sub- 
plots had higher mean ANPV and sub-plot mortality in 
comparison to low disturbance and undisturbed forest (Table 1). 
Plateaus and slopes from the high disturbance also had higher 
neighboring mortality and gap size than plateaus and slopes from 
the low disturbance. The disturbed sub-plots were set in gaps that 
ranged in size from 0-22.6 ha (6.5 ±2.3), which correlated 
positively with mean sub-plot mortality (^ = 0.37; P<0.001; 



Pearson's r = 0.60). During the analyzed period we did not find 
gaps larger than three pixels (2700 m 2 ) in the undisturbed forest. 

In the disturbed forest the density of dead trees ranged from 0 to 
360 trees ha -1 and the mean (101±22 trees ha -1 ) varied among 
topographic classes (^ = 8.74; P<0.001). The mean density of 
dead trees in plateaus (104 ±43 trees ha -1 ) was similar to that for 
slopes (120±31 trees ha ) [P = 0.536), but exceeded the mean 
density observed in valleys (44 ±31 trees ha ') (P<0.01, both). 
The total number of dead trees, both snapped and uprooted, 
correlated positively with sub-plot mortality in slopes in plateaus. 
In the disturbed forest plateaus and slopes, the density of live trees 
and basal area were negatively correlated with sub-plot mortality, 
while valleys were not affected (Table 2). 

Selective mortality patterns 

Our data show that the blowdowns promoted subsequent 
changes in forest structure, genera substitution and species 
demography. Four years after the blowdown event, differences 
in size-distribution, genera importance and diversity measures still 
correlated with levels of disturbance. 

In the disturbed forest, the DBH of recorded dead trees ranged 
from 10 to 120 cm, with mean value (27.5 ±2.3 cm) higher than 
from that from live trees of both disturbed (20.6±2.3 cm) and 
undisturbed forest (21 ±0.7 cm). DBH of dead trees did not vary 
among modes of death (ANOVA F=2.37; P = 0.125), but for 
snapped and uprooted trees had a low correlation with sub-plot 
mortality (r 2 = 0.02; P = 0.029; Pearson's r = 0.13), with a non- 
variable mean value among modes of death (ANOVA F = 2.37; 
P = 0.125). 

Even four years after the blowdown event, the DBH distribution 
comparisons indicated size distribution differences between trees 
of the disturbed and the undisturbed forests. The DBH 
distribution of dead trees in both low- and high-disturbance sub- 
plots (Figure 2a) did not follow that of the undisturbed forest 
before (2004) and after (2006) the blowdown event (Chi-squared 
test, P<0.001). While the undisturbed forest had a higher 
abundance of dead trees within the smaller diameter classes, low 
and high disturbance forests had higher abundance of large trees 
exhibiting a log-normal distribution. The DBH distribution under 
low and high disturbance was also different (Chi-squared test, 
P = 0.017). Within live trees, both disturbed and undisturbed 
forests had higher abundance of small trees and exhibit the typical 
negative exponential DBH distribution of tropical forests (Fig- 
ure 2b). The DBH distribution of live trees under low disturbance 
was similar to that of the undisturbed forest (P = 0.490), while 
under high disturbance forest it diverged from that observed 
within trees from the undisturbed and low disturbance forest (P< 
0.001). The overall mean DBH of live trees was negatively related 
to sub-plot mortality in plateaus and slopes of the disturbed forest 
(Figure 2c). 

In the disturbed forest, recorded trees were distributed in 51 
families, 158 genera and 324 species or morphotypes. A total of 54 
species (~ 16%) just occurred on plateaus, while 55 (~17%) and 20 
(~6%) were restricted to slopes and valleys, respectively. Fifty 
eight species (~ 18%) occurred in all topographic classes. Fabaceae 
(66 species), Sapotaceae (26) and Lecythidaceae (23) were the 
richer families in number of species and together summed ~40% 
of all recorded trees and ~36% of the total richness. Eighteen 
families were represented by just one species. 129 species were 
represented by only a single tree. Heavily damaged sub-plots, 
especially those in gaps with large amounts of necromass, often 
had high density of lianas. Field observations showed that gaps in 
which wood decay was more advanced generally had a dense 
understory with seedlings dominated by the genera Casearia and 



PLOS ONE | www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | e103711 



Wind Disturbances Promote Tree Diversity in Central Amazon 













o 






o 


d 




ro 


+1 


+ 1 


+i 


+ 1 






vO 






d 


d 





rsj 


00 


LO 


d 


rM 


fN 


+ 1 


+1 


+ 1 




O) 


CN 


CO 


LO 


ro 



00 


O) 


ro 










q 


d 






d 


d 


d 




fN 


+1 


+ 1 


+i 


+ 1 


+1 


+ 1 


+i 


+ 1 






rM 


q 


<*■ 


CN 


00 


CN 


ro' 






CM 




LO 




m" 



i— m rM 



rsi 


m 




o 


ro 


rsi 


o 


o 


o 


d 


d 


d 


+ 1 


+1 


+ 1 


CO 




vo 


-sT 


LT) 


-sT 


O 


o 


O 


d 


d 


d 



















Oi 




^J- 


fN 






rsi 


ro 


CN 


d 


+ 1 


+ ■ 


+ 1 


+ 1 


+1 


+ 1 


r-- 




00 


rM 


IV 


o 


r< 




<vD 


CO 


vO 


CO 




rsi 


CM 


rM 


fN 


fN 



Q\ ro rjv vD CO 



«- tTi 



,- ct» «- 











fN 


ro 


00 


CO 




LO 


VO 




+ 1 


+ 1 


+1 


+1 


q 


o 


r». 


CO 




■j5 




vd 


fN 


rsi 


(N 













fN 


r-s 


vq 






vd 


in 




+ 1 


+ 1 


+1 


+ 1 


CT> 




un 




CO 


d 


00 


CO 


fN 


ro 


fN 





ro ro ro rM 













ro 






fN 


o 


lo 




CO 




rv 


+ 1 


+ 1 


+1 


+ 1 


+1 


ro 




LO 


un 


LO 


LO 


oS 


LO 






rsi 




fN 




rM 



CO <N 00 vO U"i 

vo ov ^- 
+1 +1 +1 +1 +1 



ro = ™ 



Q 
X 



o eg I 

E.|gS 

o ° V u 
"•Da. u 

# JZ — ■ = 
m C7i vu q. 

3 f | & 

avi^ o 
3ZJj 

u ^ 0J 
c S, £ 

|j° S 

Ills 

■— ^ CL> <u 
^3 O _ g 

° s § -° - 

9 i, "= 
«■ 3 < 



S 2 

LO <— 



£ LO 



■ C oj 

<u o 

u Q_ 

c — 

<U ro 



PLOS ONE | www.plosone.org 



6 



August 2014 | Volume 9 | Issue 8 | e1 0371 1 



Wind Disturbances Promote Tree Diversity in Central Amazon 



0> !~M vO 0> 0> m 



O ,— 

V d 



i — ro lt) 

O m m 
O *- i- 



in in lo co 

O <- 



CO vD CO CO 



0 

z 



o> rsi <— 



CO 

r-. 




in 
o 


m 
o 


o 


s 


d 


d 


d 




LT) 


un 




q 


in 
o 
m 


d 




d 






ON 


m 


01 





00 LT) ^ 



01 
0 

a 
o 



S a. «i a a o) g 



Q. & Si 
° "5 JS 



Q. ^ 
° "5 

ul > 



e q- 

o 73 



" OJ 



PLOS ONE | www.plosone.org 



7 



August 2014 | Volume 9 | Issue 8 | e103711 



Wind Disturbances Promote Tree Diversity in Central Amazon 



0.4 -| 



0.3 - 



~ 0.2 - 



0.1 



0.0 - 1 



12.5 



17.5 



22.5 



27.5 



32.5 



37.5 



■ high disturbance 

■ low disturbance 
undisturbed 2004 

□ undisturbed 2006 



42.5 



47.5 



(a) 



>50 



Diameter classes (cm) - dead trees 



c 



ro 
2 



0.4 
0.3 
0.2 - 
0.1 - 
0.0 - 




■ high disturbance 

■ low disturbance 
□ undisturbed 2008 



( b ) 



x 

CD 
Q 



•a-plateau 
© slope 




12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 52.5 57.5 >60 



0 10 20 30 40 50 60 70 



Diameter classes (cm) - live trees 



Sub-plot mortality (%) 



Figure 2. Tree size distribution (DBH>10 cm) in plateaus and slopes of a forest disturbed by a large blowdown and an undisturbed 
forest, in Amazonas, Brazil. Legend: (a) diameter distribution of dead trees; (b) diameter distribution of live trees; (c) mean DBH of live trees in 
disturbed plateaus (^ = 0.05; P = 0.06; Pearson's r= -0.277) and slopes (r 2 = 0A9; P<0.001; Pearson's r= -0.451) related to sub-plot mortality. 
doi:1 0.1 371 /journal.pone.01 0371 1 .g002 



Laetia (Salicaceae), Cecropia and Pourouma (Urticaceae), Con- 
ceveiba and Croton (Euphorbiaceae), Inga and Tachigali (Faba- 
ceae), Miconia (Melastomataceae), Tapirira (Anacardiaceae) and 
Vismia (Hypericaceae). Genera demographic analyses revealed 
differences between low (total of 258 species) and high disturbance 
(215). Ten genera among the 20 most important genera from the 
high and the low disturbance levels were different (Figure 3). The 
high disturbance level had both typical light-demanding/fast- 
growing (e.g. Cecropia, Inga and Pourouma) and shade-tolerant/ 
slower-growing (e.g. Brosimum, Pouteria, Sloanea and Vantanea) 
exclusive genera within the 20 most important genera. Eschwei- 
lera, Licania, Pouteria and Protium had the highest IVI in both 
disturbance levels, but the contributions of the IVI estimators of 
these and other common genera varied between disturbance 
levels. 

Species richness per sub-sub-plot varied from 3 to 25 species in 
the disturbed forest. For groups of plots, species richness was 
negatively related to sub-plot mortality (r 2 = 0.39; P<0.001; 
Pearson's r = —0.64) (Figure 4a). Shannon diversity also correlated 
negatively with sub-plot mortality (r 2 = 0.39; P = 0.002; Pearson's 
r= —0.61). Nonetheless, at intermediate disturbance (20-50% of 
tree mortality), Shannon diversity was slightly lower or even 
similar to that from the less- and non-disturbed sub-plots 
(Figure 4b). Species rarefaction curve along the mortality gradient 
indicated that intermediate-disturbance areas can be more diverse 
than both heavily and undisturbed areas (Figure 4c). 

Disturbance gradient effects and vegetation responses 

Of the live trees recorded in the disturbed forest, 206 (~1 1% of 
total sampled) had crown or trunk injuries and 191 (~10%) had 
one of the mechanisms indicative of resprouting. 151 trees (~7%) 
belonged to light-demanding/fast-growing pioneer genera, and 
the mortality gradient appeared to amplify forest niches that lead 
to changes in species composition. 



Resprouting trees were recorded in all the topographic classes 
and within different species. In slopes, the fraction of resprouters 
related positively to sub-plot mortality as well as to other 
disturbance intensity measures (Figure 5). In plateaus and slopes, 
the mean wood density correlated negatively with sub-plot 
mortality (r 2 = 0.03; P<0.035; Pearson's r= -0.19), but resprou- 
ters and non-resprouters had no significant differences in wood 
density {F = 0.262, P = 0.609), but did have in DBH (F = 4.87, 
P = 0.027). 

In general, light-demanding/ fast-growing genera (e.g. Cecropia, 
Conceveiba, Inga, Casearia, Miconia, Pourouma, Tachigali, 
Tapirira and Vismia) contributed less to total basal area in less 
disturbed sub-plots, while their basal area was positively related to 
sub-plot mortality in plateaus (7^ = 0.39; P<0.001; Pearson's 
r = 0.63) and slopes (^ = 0.27; P<0.001; Pearson's r = 0.52) 
(Figure 6a). In contrast, the fraction of basal area contributed by 
these fast-growing species was negatively related to sub-plot mean 
wood density, also exclusively in plateaus (r 2 = 0.30; P<0.001; 
Pearson's r=— 0.56) and slopes (r 2 = 0.31; P<0.001; Pearson's 
r=— 0.57) (Figure 6b). The fraction of these typical pioneer 
species was positively related to neighboring mortality in plateaus 
(r 2 = 0.32; P<0.001; Pearson's r = 0.58) and slopes (r 2 = 0.46; P< 
0.001; Pearson's r = 0.67) (Figure 6c). Gap size was also positively 
related to these group of species in plateaus (r 2 = 0.14; P = 0.013; 
Pearson's r = 0.34) and slopes (r 2 = 0.21; P<0.001; Pearson's 
r = 0.44) (Figure 6d). 

The NMDS stress value was 0.286, and the first two axes 
captured 34% of the floristic variation (Figure 7a and Table 3). 
Floristic similarities within heavily damaged sub-plots increased 
due to the higher abundance of light-demanding and fast-growing 
species. Sub-plot mortality related positively to the variations in 
species composition captured by the first NMDS axis (Figure 7b). 
The same pattern was observed for neighboring mortality and gap 
size, which correlated more strongly (Table 3). 
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Figure 3. Genera importance ranking over plateaus and slopes of a forest disturbed by a large blowdown, in Amazonas, Brazil. 

Legend: (a) Importance Value Index (IVI) of the 20 most important genera recorded under low disturbance (up to 13% of tree mortality); (b) IVI of the 
20 most important genera recorded under high disturbance (up to 70% of tree mortality). 
doi:1 0.1 371 /journal.pone.01 0371 1 .g003 



The demographic patterns of the 25 most important genera in 
the community (identified in Figure 3) and the variations in species 
composition (Figure 7) corroborate our results related to genera 
abundance variation in respect to disturbance intensity (Figure 8). 
The genera abundance curves show that these 25 genera are not 
equally distributed along the landscape and that part of this 
variation was regulated by variations in sub-plot mortality. 
Interestingly, the abundance of these genera varied along the 
mortality gradient revealing five distinct groups (identified using 
colors in Figure 8) which had their optimum under specific 
disturbance intensities. This pattern indicates that there exist 
specialist guilds for each different level of disturbance/ resource. 
Moreover, areas where the tree mortality driven by the blowdown 



ranged between 20-50% are those which provide gap space and 
light conditions, thus favoring a wider range of species with 
different requirements. 

Discussion 

Disturbance patterns 

Mortality rate and structural measurements of the undisturbed 
forest were similar to those reported in other undisturbed forests 
nearby [52,67,77,85]. Tree mortality levels observed in the 
disturbed forest were in accordance with patterns described in 
earlier studies [9,10,25] also in Amazon and even higher than in 
hurricane-damaged forests [30,45,47,86]. 




Figure 4. Species richness and diversity measures related to a mortality gradient in plateaus and slopes of a forest disturbed by a 
large blowdown, in Amazonas, Brazil. Legend: (a) species richness; (b) Shannon diversity; (c) species rarefaction. 
doi:1 0.1 371 /journal.pone.01 0371 1 .g004 
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Figure 5. Fraction of resprouting tree species related to disturbance intensity measures from a forest disturbed by a large 
blowdown, in Amazonas, Brazil. Legend: (a) fraction of resprouters related to sub-plot mortality (slopes Pearson's r = 0.33); (b) fraction of 
resprouters related to neighboring mortality (slopes Pearson's r = 0.38); (c) fraction of resprouters related to gap size (slopes Pearson's r = 0.33). 
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Wind-related mortality caused by snapping and uprooting was 
greatest in more exposed areas, including plateaus and the top of 
slopes (Table 2). As a consequence of lower disturbance intensity, 
structural patterns were not altered by blowdowns in the valleys 
(Table 1). This may indicate that strong winds did not reach the 
valley floors, or that the vegetation in valleys is less affected by 
wind damage, as has been observed in other tropical and sub- 
tropical forests affected by cyclones [23,24,29,48,87]. We also did 
not observe uprooted canopy palm trees in valleys and just 14 
snapped and/or standing dead ones. As there is no evidence that 
caulescent palm tree species that occur in these valleys are less 
susceptible to wind damage (see species in item 2.2), we assume 
that valley floors are indeed better protected from strong winds. 



Although in these forests vegetation damage was partially 
controlled by abiotic aspects such as wind characteristics and 
topography [22,29,88], resolution of this question requires more 
information on individual-species traits, such as anatomic- 
morphological variations among individuals and species, size 
distribution and populations range, and individual pre-disturbance 
conditions (pest-attack, biological interactions, age, etc.). 

In the disturbed forest, measures of forest structure, including 
tree density and basal area, were lower four years following the 
blowdown event (Table 1). The topographic variation of measures 
of structure was amplified by the landscape driven differences in 
mortality. Immediate tree mortality effects on tree density and 
basal area reduction are direcdy dependent on the disturbance 




Figure 6. Pioneer (light-demanding and fast-growing) tree species importance related to disturbance intensity measures and wood 
density in a forest disturbed by a large blowdown, in Amazonas, Brazil. Legend: (a) fraction of pioneers' basal area related to sub-plot 
mortality; (b) fraction of pioneers' basal area related to sub-plot mean wood density; (c) fraction of pioneers related to neighboring mortality; (d) 
fraction of pioneers related to gap size. 
doi:1 0.1 371 /journal.pone.01 0371 1 .g006 
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Figure 7. NMDS ordination diagram of 1 17 sub-plots sampled in plateaus and slopes of a forest disturbed by a large blowdown, in 
Amazonas, Brazil. Sub-plot mortality caused by the blowdown is highlighted by the color-scale. Legend: (a) NMDS scores computed from a 
dissimilarity matrix with the abundance of all recorded species; (b) first NMDS axis related to sub-plot mortality. 
doi:1 0.1 371 /journal.pone.01 0371 1 .g007 



intensity - i.e. the fraction of trees killed by wind (Table 2). In 
heavily damaged areas with the highest estimated mortality rates, 
structural measures were lower than the early-successional stages 
found in other nearby and smaller disturbed forest patches [89] 
and anthropogenic disturbances in the same forest type [41,90]. At 
these sites, forest recovery will depend on species responses to the 
light gradient in gaps, which was supported for our results related 
to the third and fourth questions. Within the disturbed areas, the 
correlation of percent mortality with structural characteristics also 
showed that areas less affected by wind damage were similar to the 
undisturbed forest sub-plots. 

The observed higher mortality within plateaus and slopes of the 
undisturbed forest was also related in nearby unaffected terra- 
firme forest, where under regular disturbance regimes, standing 
death can be expected to be more frequent in plateaus, and 
uprooting and snapping more frequent in slopes and in valleys 
[91]. Thus, competition can be more important in plateaus and 
exogenous disturbances dominant in slopes and valleys, which will 
cause differences in the size-density distribution of trees related to 
topographic position [92]. Despite such findings, our study 
indicates that for the Central Amazon forests, most notably in 
plateaus and slopes, tree mortality is also regulated by wind 
disturbance regimes such as blowdowns of varying size. 

Although in the undisturbed forest we only observed small gaps, 
in the disturbed forest large gaps ( 2000 m 2 ) were common and 
changed forest structure, particularly in plateaus and slopes. Such 
results suggest that severe disturbances are required to form large 
gaps. The higher values of tree density and basal area observed in 
plateaus may indicate a variation in forest vulnerability due 
differences in species composition, or possibly higher resilience of 
plateaus. 

Selective mortality patterns 

We expected to see large emergent to be more prone to wind 
damage than smaller ones. While we did observe differences in 
DBH distribution of dead trees between undisturbed and 
disturbed forests that suggests lower proportional mortality for 
smaller trees (despite potential biases introduced because smaller 
dead trees decompose faster), we also observed no correlation 
between the DBH of dead trees and sub-plot mortality, indicating 



that size selective mortality did not depend on disturbance 
intensity. 

Observed reductions in mean DBH and wood density in the live 
trees in highly damaged areas of the disturbed forest (Figure 2c) 
likely reflect fast recruitment of light-demanding and fast-growing 
species with shorter life span and generally lower wood density 
values [15,34,71]. Chambers et al. (2009b) compared non- 
disturbed and smaller blowdown gaps with ~10 years-succession 
and reported similar patterns, except for tree density, which was 
lower in our four-years disturbed plots. The substantial increase in 
the fraction of small trees with increasing sub-plot mortality also 
evidentiate the recruitment of shade-tolerant or understory species 
that usually do not reach the canopy, and thus demonstrates a 
selective mortality feedback on floristic composition. 

Tree diversity in Central Amazon forests is high, and our data 
agree with those of other studies conducted in the same [52,67,93] 
and in adjacent areas [55,57,94,95]. We observed clear effects of 
the blowdown mortality on community demographic patterns 
from low to high disturbance (Figure 3). Part of the observed 
changes can be attributed to the fast recruitment of species with 
different syndromes and growth (lateral expansion) in the available 
gap space, such as Cecropia sciadophylla Mart., Pourouma 
tomentosa Mart, ex Miq., Inga cf. paracusis, Iryanthera juruensis 
Warb., Protium hebetatum D. C. Daly, Mabea speciosa MtilLArg. 
and Zygia racemosa (Ducke) Barneby & J. W. Grimes. These 
species improved their importance within the high disturbance 
sub-plots by increasing their relative density and frequency. Such 
patterns may indicate success of establishment and expansion of 
these populations [81], which in the longer-term may increase 
species richness and forest diversity [89]. Common tree canopy 
genera in adjacent areas (e.g. Couratari, Eschweilera, Pouteria, 
Licania, Micropholis, Ocotea, Sloanea, Swartzia and Vantanea) 
[67,90,93], also presented high importance within the sub-plots 
with high disturbance (i.e. high mortality). Survivorship of shade- 
tolerant and slower-growing species from these genera may be 
attributable to the relatively high wood density (~0.79 g/cm 3 ) 
[71], small crown footprint, and presence of buttress and 
supporting roots [54,56]. Blowdowns are a major disturbance in 
these forests, thus the greater likelihood of survival and resprouting 
ability may be reasons for these genera to be of such high 
importance in the terra firme forests of this region [54—56,96] . 
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Blowdown mortality, which direcdy affected species richness by 
reducing tree density at a local scale, resulted in the negative 
correlations between species richness and Shannon diversity with 
sub-plot mortality (Figure 4a and b). Nevertheless, the rarefaction 
curve indicate that a slightly higher number of species would be 
observed for any smaller subsample of individuals taken from low 
to intermediate disturbance areas, under the assumption of 
random mixing of individuals (Figure 4c). Interestingly, despite 
differences in species richness driven by disturbance intensity, the 
higher number of rare species at low and intermediate disturbance 
levels indicates that a short-term feedback of blowdown mortality 
is to promote species richness at the landscape scale, with apparent 
lower impacts on valleys. 

Disturbance gradient effects and vegetation responses 

So far we have focused on the importance of wind disturbance 
for providing niche space for light-demanding species requiring 
large canopy gaps for regeneration as a mechanism for promoting 
species coexistence and thus diversity at the landscape scale. 
Concomitantly, our data and tests related to the third and fourth 
questions support that the studied gaps also provide adequate 
niches for species with different requirements. 

Disturbance and selective mortality patterns promoted imme- 
diate effects on forest structure and species composition, as 
evidenced by changes in tree size distribution (Figure 2), species 
importance between levels of disturbance (Figure 3) and species 
diversity along the mortality gradient (Figure 4). Nonetheless, 
blowdown gaps seem to produce long-term effects at the 
community level. The observed landscape mortality gradient also 
changed diversity patterns by promoting compositional changes 
with respect to the fraction of resprouters (Figure 5) and with light- 
demanding/ fast-growing species (Figure 6). 

Thus, a second potentially important niche axis is the ability to 
survive wind disturbance, which may provide a selective advan- 
tage to species that are less susceptible to wind damage either 
because they possess specific morphological and anatomical 
adaptations lending stability or because they are of small stature 
and are typically not exposed to wind. Both characteristics may 
g. represent a competitive disadvantage in the absence of wind 

disturbances. In the first case, stability conferred by e.g. high wood 
density is associated with high construction costs and reduced 
growth rates - a typical growth-defense tradeoff reported in the 
literature [8,71,97]. In the second case, staying short or retaining a 
low height/diameter ratio is certainly a disadvantage in the race 
for light [62]. 

The positive correlations between disturbance intensity mea- 
sures and density of resprouters in the disturbed forest (Figure 5) 
confirms the importance of this regeneration mechanism in gaps 
created by wind disturbances [31,47,48,86,98,99]. The similarities 
in mean wood density between resprouters and non-resprouters 
(usually undamaged trees) indicate that resprouting is a regener- 
ation pathway adopted by species with different traits and light 
requirements, predominandy by smaller trees (DBHS20 cm). A 
similar pattern was found in a Caribbean hurricane-damaged 
| <s\ jjj forest [45] and indicates that, although canopy trees from higher 

wood density species may have higher wind resistance, resprouting 
is not an exclusive regeneration pathway for shade-tolerant, 
slower-growing or climax species. As for sub-tropical hurricane- 
damaged forests [46] , we hypothesize that in terra firme forests of 
Central Amazon, both early secondary and understory species 
may benefit from wind disturbances. Finally, the decrease in the 
mean aggregated wood density in disturbed sub-plots is a direct 
consequence of the higher fraction of pioneer species which are 
characterized by low wood densities (Figure 6). Wood density 
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Sub-plot mortality (%) 

Figure 8. Genera abundance related to a mortality gradient in plateaus and slopes of a forest disturbed by a large blowdown, in 
Amazonas, Brazil. Legend: low disturbance specialist genera- Brosimum, Eperua, Iryanthera, Micrandropsis, Micropholis, Pouteria, Protium and 
Unonopsis; resistant or unaffected genera- Couratari, Ocotea and Sloanea; genera favored under intermediate disturbance- Eschweilera, Lecythis, 
Licania, Naudeopsis, Swartzia, Qualea and Zygia; genera favored under high disturbance- Couepia, Mabea, Sderonema and Vantanea; high disturbance 
specialist genera- Cecropia, Inga and Pourouma. 
doi:1 0.1 371 /journal.pone.01 0371 1 .g008 



variation among sub-plots was partially explained by sub-plot 
mortality, which is significant when low wood density pioneers are 
observed (Figure 6a and b). This greater importance of pioneers in 
heavily damaged areas corroborates demographic patterns found 
by assessing species IVI (Figure 3b) and patterns in species 
composition, which was indicated by the NMDS-test and 
regressions (Figure 7). 

Our data show that species turnover in the 20 most important 
genera (as quantified by IVI in Figure 3) is due to the increasing 
importance of both non-pioneer and pioneer species. Indirect 
evidence for this mechanism is provided by fact that the floristic 
dissimilarity between undisturbed and disturbed sub-plots seems to 
be driven by the admixture of light-demanding species (Figure 7a 
and Figure 8). There was a significant correlation between sub- 
plot mortality and the position of sub-plots along the first NMDS 
axis, the latter explaining significant part of the total variation in 
species composition (Figure 7b). This relationship was higher 
when considering neighboring-mortality and gap size values 
(Table 3). However, even in the plots with highest mortality rates 
(40-70%) there was no complete takeover by classical pioneer 
local genera (e.g. Cecropia, Conceveiba, Croton, Goupia, Inga, 
Laetia, Miconia, Pourouma, Tachigali and Vismia) or shade- 
tolerant common genera in undisturbed forests (e.g. Couratari, 
Eschweilera, Licania, Pouteria, Sderonema, Sloanea and Swart- 
zia). These last, maintained up a large fraction of the IVI as 
indicated by their high proportion in the largest diameter class. 
Some of the fast growth species we have recorded, especially from 
the genus Cecropia, Inga, Pourouma and Tachigali were also 
reported as important species during early succession stages of 
secondary terra-firme forests around Manaus [40,90,100] and in 
the upper Rio Negro [41]. In anthropic secondary forests, under 
non-intensive use and without fire regimes, areas dominated by 
species from the genus Cecropia and Pourouma present relatively 
higher species richness, which hence more rapid plant succession 
[40,100]. 

The variations of genera importance in respect to mortality 
intensity support classical studies [2,5,6,34] that attempt to classify 
species into a simplified conception of pioneer species in gaps 



versus intermediate or late-succession species under undisturbed 
forest patches. These classical opposing strategies were observed 
within the 25 most important genera of the disturbed forest 
(Figure 8), as characterized as high (light-demanding and faster- 
growing - red lines) and low disturbance specialists (shade-tolerant 
and slower-growing - green lines). Additionally, our data 
highlighted alternative successional trajectories related to survival 
(resistant or unaffected - black lines), resprouting and fast- 
recruitment (favored under intermediate and high disturbance - 
blue and orange, respectively). 

These alternative successional trajectories have not been 
observed as important regeneration mechanisms in smaller treefall 
events. In larger gaps (>2000 m 2 ), the canopy emergence of 
specialized guilds with an optimum at specific levels of 
disturbance, indicate a smooth occupation of the entire gap space 
and light conditions. The similar pattern within genera belonging 
to the same guilds may indicate similar regeneration strategies and 
historical life. A possible explanation for the success of species from 
these different genera may be the higher plasticity of traits and 
alternative resilience mechanisms (such as survival, resprouting 
and recruitment). As blowdowns are common in these forests 
[10,12,14], these results show that a simplified classification system 
is not enough to describe successional trajectories of large gaps. 

Surprisingly, all the 25 most important genera in the disturbed 
forest (Figure 3) also figured within the 1 2 1 most abundant genera 
of the main portion of the Amazon [96] . Protium, Pouteria and 
Eperua characterized the low disturbance specialist guild in the 
disturbed forest (Figure 8) and took the second, third and seventh 
position in the Amazonian genera abundance ranking, respectively 
[96]. Ocotea, Sloanea and Couratari characterized the resistant or 
unaffected guild in the disturbed forest and took the 14°, 22° and 
106° position in the Amazonian ranking, respectively. Eschweilera, 
Licania, Lecythis and Swartzia characterized the intermediate 
disturbance intensity at the disturbed forest and took the first, 
fourth, 10° and 17° position in the Amazonian ranking, 
respectively. Inga, Pourouma and Cecropia characterized the 
high disturbance specialist guild in the disturbed forest and figured 
at the sixth, 18° and 30° position in the Amazonian ranking. 
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The co-occurrence of species with a broad range of life history 
strategies at intermediate disturbance levels (evidentiated as higher 
confluence of lines - Figure 8) shows that this blowdown promoted 
species richness probably by amplifying niches and/or resources. 
We assume that large gaps (>2000 m 2 ) produced by blowdowns 
can be pardy associated to increased diversity and the contrasting 
dominance of some species in these forests [54,55,57,67,94,96]. 
Variations in genera importance and species composition along 
the disturbance gradient indicate that gap recovery in large gaps is 
influenced by mortality patters, most pronounced in plateaus and 
slopes. Additionally, selective mortality and fast-recruitment 
highlight genera specialization along the mortality gradient, and 
indicate that large gaps contain environmental variability that 
together with species responses and resistance, allow high tree 
species diversity. In this study, intermediate-disturbance levels had 
higher species richness suggesting that there was a selective 
mortality and resulting species turnover in response to the 
disturbance gradient. This pattern provides evidence that blow- 
downs allow species richness and diversity through an interaction 
of wind-damage with species-resistance and resilience. Although 
classical successional guilds may dominate small treefall gaps, the 
observed pattern indicate that in large gaps (>2000 m 2 ) there are 
specialists for a wider range of disturbance and light. Thus, we 
assume that a diverse set of species differing widely in light 
requirements and recruitment strategies forms the pioneer cohort, 
thus lending a high resilience towards wind disturbances at the 
community level. 

Our results add new and complementary information about 
succession and turnover of Western and Central Amazon forests 
and reinforce that recovery processes for large gaps differ from the 
trajectory observed in gaps formed by smaller treefall events 
[2,6,15,31,86] and secondary forests from anthropic activities, 
where fire and logging may increase biomass losses [41,90] and 
limit species regeneration [40,100]. Moreover, we revealed a 
gradient in demographic responses along the entire gradient of 
disturbance levels. Mortality intensity and gap size seem to 
influence community composition by filtering shade-tolerant 
survivors with resprouting ability and favoring more light- 
demanding and fast-growing species. As already observed in 
tropical and subtropical regions [26,29,32,47,98,101] wind dam- 
age depends on species composition and successional stage, which 
suggest that secondary forests in the Amazon might be more 
vulnerable and less resilient to windstorms than forests in more 
advanced successional stage. Thus, for better understanding the 
blowdown effects on species distribution, forest vulnerability, 
ecosystem functioning, and the ecological importance of large gaps 
on species maintenance it is necessary that future studies include 
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